Microscopic Density Functional Theory for Dendrimers 
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Density functional theory for a simple model of dendrimers is proposed. The theory is based on 
fundamental measure theory which accounts for the hard-sphere repulsion of the segments and on 
the Wertheim first-order perturbation theory for the correlations due to connectivity. Set of the 
recurrence formulae for the ideal chain contribution involving simple integrals is derived. By using 
perturbation theory dispersion forces can be easily included. 



Dendrimers (also known as arborols or cascade 
molecules) are repeatedly branched monodisperse com- 
pounds with a fractal-like structure possessing a high 
degree of symmetry. The first dendritic structure was 
constructed by Vogtle et al. in 1978 [l[ using a repet- 
itive synthesis strategy (divergent synthesis) such that 
one new molecular layer (generation) is created in each 
reaction cycle. Almost immediately dendrimers attracted 
increasing attention due to their unique structure imply- 
ing some unique properties compared to traditional lin- 
ear chains. With increasing number of generations the 
number of chain-ends increases exponentially and that is 
why dendrimers adopt a compact globular shape. As a 
consequence the dendrimers' solubility is driven by the 
nature of surface groups only Moreover, the pres- 
ence of internal cavities enables encapsulation of small 
guest molecules Q. Such properties predestinate den- 
drimers for a variety of possible technological application. 
They are being developed for use in fields such as cataly- 
sis, magnetic resonance imagingjdrug delivery, coatings, 
electronics or cancer therapy [4j-|7[. 

In terms of theoretical study and in particular of phys- 
ical model used, there are in principle two perspectives 
in treating such complex structures. First, one can use 
so-called coarse graining technique treating molecules as 
spherically symmetric objects interacting via some ef- 
fective soft potential, i.e. the degrees of freedom of 
monomers composing the polymer are integrated out. 
Perhaps the first who proposed such a strategy was F lory 
8], but it has become popular only in recent times jj^lioj. 
The coarse-grained approach was applied for dendrimers 



by Gotze et. al [ll| who approximated the interac- 



tion between two dendrimers in solution by an appropri- 
ate Gaussian function. Even though the coarse-graining 
method is attractive in simplifying the description of a 
given complex structure substantially, it inevitably loses 
some information about the intrinsic property of the 
molecules. Moreover, the radial symmetry of an effec- 
tive potential is more justified in a bulk phase rather 
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than e.g. in the vicinity of a wall. For these reasons, 
a second approach, treating the complex molecular sys- 
tems on an atomistic level, would seem to be to superior. 
Microscopic density functional (DF) theory provides a 
versatile and powerful tool to represent the microscopic 
structures and interfacial phenomena of polyatomic flu- 
ids under a variety of situations. Woodward developed 
a theory that combines weighted density approximation, 
known from theories of simp le fluids, with single-chain 
Monte Carlo simulations [12| . An alternative DFT of in- 
homogeneo ns p olymer solutions was formulated by Fors- 
man et al. [13[ . Their theory is based on the free energy 
functional resulting from the generalized Flory equation 
of state. However, an approach due to Yu and Wu [H| 
which incorporates Wertheim's perturbation theory for a 
bulk fluid [15l. Il6| into the non-local DF framework pro- 
posed by Rosenfeld [lj} revealed to be particularly ap- 
pealing. In the spirit of the Rosenfeld fundamental mea- 
sure theory (FMT), Yu and Wu constructed a non-local 
functional accounting for the chain connectivity. Such 
an approach has proved to be both quantitatively accu- 
rate and computationally convenient and eventually has 
been extended for a variety of models of e.g. cyclic poly- 
atomic fluids [lH, block copolymers [l9| . star polymers 
[jo| . polydisperse polymers [Il| or for brush-like struc- 
tures [22j . Recently, a so-called hybrid approach for the 
structure of dendrimers has been proposed [23| in the 
spirit of Ref. [12]. In this paper, a full density func- 
tional is derived and used to represent hard dendrimers 
confined between two hard walls. 

In the following, by the term "dendrimer" will be taken 
to mean the special case of a tree structure where each 
segment apart from the terminating ones has the same 
number of bonds, three at minimum. Due to the high 
level of symmetry the system can be characterized by 
two parameters: i) /, number of arms, i.e. bonds out- 
coming from each (except the terminating) segment; ii) 
M , the number of generations, i.e. number of segments 
contained in a chain connecting the central and termi- 
nating segment minus one. Dendrimers with / arms and 
M generations will be abbreviated by D(f, M). 

Each dendrimer contains a central a segment, which is 
by definition segment of generation 0. Segments of the 
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ith generation are connected to the central segment by 
a chain of i + 1 segments. The number of segments of 
ith generation is <?j = / • (/ — x , i > 1 and the total 

number of segments is N = f 
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Segment positions are labeled by two indexes; the sub- 
script, i = 0..M, specifying the generation, and the su- 
perscript, j = l..g(i), specifying the position in a given 
generation. The latter can be set in a clockwise order 
(in a two-dimensional projection) such that segments 
r rj" 1 are connected to the segment r^. Posi- 
tion of a whole dendrimer can be expressed by a vector 

The model under interest will be represented by tan- 
gentially connected hard spheres of diameter a, each in- 
teracting via potential ip with an external field. The 
grand potential functional of such a system can be ex- 
pressed as 



pn[p N (K)} = (3F id [ PN (R)]+l3F cx + J [*(R)-/i] P7V (R)dR, 

(1) 

where 



PFj, 



dR PiV (R)[logpAr(R)-l]+/3 J dK PN (R)V b (R) 

is the contribution corresponding to the system of ideal 
chains that interact only through bounding potential, 
Vf,(R), and the excess part that takes into account cor- 
relations between nonbonded segments 



/^cx = J dr{$ hs K(r)] +$ c K(r)]} 



(3) 



is split into the hard-sphere contribution $ hs and the 
contribution due to the chain connectivity $ c . 

Piv(R) is the dendrimer density and ^(R) = 

J2ito J2j=\ ^'('"i )- Further, H(R) is a sum of bounding 
potentials between the neighboring segments creating the 
dendrimer structure, 
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(4) 

Free energy densities $ hs [n Q (r)] and $ c [n Q (r)] arc 
functions of four scalar and two vector weighted densi- 
ties {n a (r)} [l4, 17 1 . For the hard-sphere contribution, 
$ hs [n a (r)] , the so-called White-Bear approach (or mod- 



1J, [2J] for the 



ified FMT) have been used, see Refs. 
explicit formulae. 

The free energy density due to indirect chain con- 
nectivity, $ c [n Q (r)], was obtained as a generalization 
of Wertheim's first-order perturbation theory for a bulk 
fluid PHH for inhomogeneous systems within the non- 
local DF framework [lij ] 



1 — N 

$> Q (r)] = -^—n (Hg HS (a,{n a })] , (5) 

where £ = 1 — ri2 • ri2/(n2) 2 , and the contact value of 
the hard-sphere pair correlation function, <?hs(ct, {n a }), 
is obtained from the Carnahan-Starling equation of state. 
The important feature of this approach is that the prob- 
lem is formulated on the level of average segment density, 
p(r), which is related to the density of a whole dendrimer 
pjv(R) via 
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where pi (r) is the density distribution of an individual 
segment. 

A minimization of the grand potential functional with 
respect to the density distributions gives rise to a set of 
the following Euler-Lagrange equations: 



p>.(r) = exp(/fy) J dRS (r - rfj cxp [— j8V&(R)] 7 (r) 



(7) 



and 
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Due to symmetry, pj (r) depends on its generation num- 
ber only, so that the upper index will be omitted. Func- 
tion 7(1") is defined as 
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Specifically, the segment density distributions have 
been calculated for hard-sphere dendrimers confined by 
two plane hard walls placed a distance H apart, i.e. the 
external field interacts with each segment with a poten- 
tial 



■ oo z < a/2 or z > H — a/2 . 
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For this system pj(r) = p\(z), and the the Euler- 
Lagrange equations have much simpler forms: 



p (z) = exp(Pvh(z){G M (z)) 
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(11) 



and 



Pi(z) = exp(^) 7 (^)(G Af _ l (z)) / - 1 G l (z) i > 1 . (12) 
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FIG. 1: (a) The average segment density of third generation 
dendrimers, D(3, 3), and linear chains, of N = 22 segments, 
(b) Segment density profiles of the first and eleventh segments 
of linear chains, of N = 22 segments, (c) Segment density 
profiles of the zeroth, first, second, and third generations of 
the dendrimer D(3, 3). The bulk segment density is pi = 0.5. 
The hard walls are separated by H = 10a. 



The functions Gi(z) and Gi(z) are defined by the follow- 
ing recurrence relations 



Gi(z) = / dz'7(;0(Gi_i(;0) 
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and 
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with Gq(z) — 1 and 



Gx(z) = / d/ 7 (2')(G M (^)) 
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In Fig. 1 comparison of the density profiles for third 
generation D(3, 3) with the profiles of linear chain poly- 
mers comprising the same number of tangent segments, 
i.e. N = 22, is shown. The latter are obtained from 
the theory of Yu and Wu [l4|]. The calculations have 
been performed for an average bulk segment density of 
Pb = Pb^ 3 = 0.5. It is evident from the upper panel 
that the average segment densities of dendrimers and 
chains are very similar for such a density. Both profiles 
exhibit oscillation characteristic with adsorption on the 
wall. The system is thus in the regime where the struc- 
ture of the fluid is dominated by excluded volume effects. 
In such a case the specific architecture of the molecules 
plays a less important role and the system as a whole be- 
haves like a hard-sphere system. It is interesting to note 
the identical contact densities, p(a/2), which reflects an 
equality of bulk pressures of both systems according to 
the sum rule [2|| . This is because the systems are treated 
in the TPT1 approximation where only number of bonds 
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FIG. 2: Density profiles of dendrimers D(3, 1), D(3, 3), and 
D(3, 5) for H — 10a and pi = 0.4. (a) The average segment 
density; (b) Segment density of cores (segments of zeroth gen- 
eration); (c) Segment density of terminal segments. 



(not their topology), which are same for both models, are 
taken into account. 

Although the average segment densities are nearly 
identical, there are large differences in the densities of in- 
dividual segments. In Fig. lb the density profiles of the 
terminal and middle segments of chains are plotted. They 
both exhibit qualitatively similar behaviour to the aver- 
age segment density, with higher value of contact density 
for the terminal segment. This is because of a smaller 
loss of orientational entropy if the terminal rather than 
\]m Middle segment is at the contact with the wall in the 
case of chains and because of excluded volume interac- 
tions from the outer segments in the case of dendrimers. 
The density profiles of the remaining segments smoothly 
interpolate between these two curves. The behaviour ob- 
served for dendrimers is very different, see Fig. lc. In this 
case only the terminal segments are in a regime where 
excluded volume effect dominates whereas all other seg- 
ments exhibit surface depletion. Clearly, the terminal 
segments can be adsorb on to the wall more easily then 
those of lower generations. Because the number of seg- 
ments of the highest generation is more than half of the 
total number of segments the adsorption for the average 
segment density persists. In all cases a cusp in the den- 
sity profiles a distance a from the sphere-wall contact 
position which reflects harshness of the fluid-fluid and 
the fluid-wall interaction. 

The impact of the generation number on a structure 
of the fluid is examined in Fig 2. The calculations are 
carried out for D(3, 1), .0(3,3), and .0(3,5) for a bulk 
density pi = 0.4. Now significant differences for dif- 
ferent architectures are apparent in the average density 
profiles, see Fig. 2a. We observe a transition from sur- 
face adsorption for £)(3, 1) (the simplest star polymer) 
to depletion for D(3, 3) and D(3, 5). Interestingly for the 
third generation dendrimer D(3, 3) the effect of deple- 
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FIG. 3: Solvation forces between two hard walls separated 
by third generation dendrimers D(3,3) (continuous curves) 
and linear chains of N = 22 segments, (dashed curves). The 
average bulk segment densities are (a) pb = 0.1, (b) pb = 0.3, 
(c) p b = 0.4, (d) p b = 0.6 



tive model of dendrimers is proposed. A compact recur- 
sive formulae is derived involving both intra- and inter- 
molecular forces taking the form of simple integrals which 
greatly facilitates the numerical calculations. This is the 
first step in a theoretical treatment of more realistic mod- 
els of dendrimers for which the theory can be straight- 
forwardly extended by using of a perturbation technique. 
It will enable to study various interesting problems; for 
instance, one of the controversy is whether dendrimers 
of higher generations adopt a membrane-like surface [26| 
or exhibit rather homogeneous segmental density due to 
back- folding effects 27 1 . Behavior of dendrimers in con- 
centrated solutions is also of an interest. The quantita- 
tive agreement of DF theory for star polymers (the sim- 
plest case of dendrimers) derived on a similar basis [1(3] 
with Monte-Carlo data gives one confidence that the the- 
ory developed in the present work provides an accurate 
representation of dendrimeric systems. However, only 
more detailed numerical tests for corresponding models 
will reveal to what extent the proposed theory for den- 
drimers is appropriate, a matter of current research. 



tion and adsorption almost compensate each other. The 
density profiles of the zeroth generation segments (the 
cores) and the terminal segment are compared in Figs. 
2b and 2c, respectively. Whereas in all three cases the 
zeroth segments exhibit depletion, adsorption is always 
found for terminal segments. Larger differences are ap- 
parent for the central segments, particularly in the slopes 
of the profiles beyond the cusps which change from neg- 
ative to positive with increasing M, being close to zero 
for D(3,3). 

In a final analysis the solvation forces between the two 
hard walls are calculated for two models, D(3, 3) and the 
equivalent linear chains of N — 22 tangent segments. 
According to the sum rule [25j, the solvation force Fs 
per unit area A is related to the average contact den- 
sity through F s /{2Ak b T) = p(0) - Poo (0) , where Poo (0) 
is the average contact density for infinite separation. In 
Figure 3 the four regimes corresponding to bulk densi- 
ties pb = 0.1,0.3,0.4, and 0.6 that the solvation forces 
can obey for the two models under consideration are 
presented. At the lowest density, Fig. 3a, the deple- 
tion forces dominate, so that the solvation force is at- 
tractive for small separations. At larger separations the 
attraction decays monotonically to zero in the case of lin- 
ear hard-sphere chains whereas for dendrimers Fs first 
changes sign and eventually converges to zero. For in- 
termediate densities, Fig. 3b and Fig. 3c, both solvation 
force profiles exhibit a maximum which is smooth for 
Pb = 0.3 but which for the higher density of pb = 0.4 be- 
comes a cusp. At the highest density of pb = 0.6, Fig. 3d, 
the specific topology of molecules becomes irrelevant and 
both solvation force profiles have very similar oscillatory 
characteristics of hard-sphere systems. 

In this work a density functional theory for a primi- 
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